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ABSTRACT 



Context. We have performed a comprehensive multiwavelength analysis of a sample of 20 starburst galaxies that show 
S~H, a substantial population of very young massive stars, most of them classified as Wolf-Rayet (WR) galaxies. 

Aims. We have analysed optical/iVZR colours, physical and chemical properties of the ionized gas, stellar, gas and 
dust content, star- formation rate and interaction degree (among many other galaxy properties) of our galaxy sam- 
ple using multi-wavelength data. We compile 41 independent star-forming regions -with oxygen abundances between 
12+log(0/H)= 7.58 and 8.75-, of which 31 have a direct estimate of the electron temperature of the ionized gas. 
Methods. This paper, only submitted to astro-ph, compiles the most common empirical calibrations to the oxygen 
abundance, and presents the comparison between the chemical abundances derived in these galaxies using the direct 
method with those obtained through empirical calibrations, as it is published in Lopez-Sanchez & Esteban (2010b). 
Results. We find that (i) the Pilyugin method (Pilyugin 2001a,b; Pilyugin & Thuan 2005), which considers the R23 and 
the P parameters, is the best suited empirical calibration for these star-forming galaxies, (ii) the relations between the 
oxygen abundance and the N2 or the O3N2 parameters provided by Pettini & Pagel (2004) give acceptable results for 
objects with 12+log(O/H)>8.0, and (iii) the results provided by empirical calibrations based on photoionization models 
(McGaugh, 1991; Kewley & Dopita, 2002; Kobulnicky & Kewley, 2004) are systematically 0.2 - 0.3 dex higher than the 
5^ ■ values derived from the direct method. These differences are of the same order that the abundance discrepancy found 

between recombination and collisionally excited lines. This may suggest the existence of temperature fluctuations in 
I the ionized gas, as exists in Galactic and other extragalactic H n regions. 

Conclusions. All these results are included in the paper Massive Star Formation in Wolf-Rayet galaxies IV. Colours, 
chemical-composition analysis and metallicity-luminosity relations, Lopez-Sanchez & Esteban (2010b), A&A, in press 
(Sect. 4.4 and Appendix A). Please, if this information is used, reference that paper and NOT this document, which 
have been only submitted to astro-ph to emphasize these results. 
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Furthermore, today it is the metallicity (which reflects the gas reprocessed by stars and any exchange of gas between 
the galaxy and its environment) and not the stellar mass (which reflects the amount of gas locked up into stars) of 
a galaxy the main problem to get a proper metallicity-luminosity relation, so that different methods involving direct 
estimates of the oxygen abundance, empirical calibrations using bright emission-line ratios or theoretical methods based 
on photoionization models yield very different values (i.e., Yin et al. 2007; Kewley & Elisson, 2008). 

Hence precise photometric and spectroscopic data, including a detailed analysis of each particular galaxy that al- 
lows conclusions about its nature, are crucial to address these issues. We performed such a detailed photometric and 
spectroscopic study in a sample of strong star-forming galaxies, many of them previously classified as dwarf galaxies. 
The majority of these objects are Wolf-Rayet (WR) galaxies, a very inhomogeneous class of star-forming objects which 
share at le ast an ongoing or recent star form ation event that has produced stars sufficiently massive to evolve into the 
WR stage (jSchaerer. Contini fc Pindad lT999V The main aim of our study of the formation of massive stars in starburst 
galaxies and the role that the interact ions with or between dwarf galax ies and/or low surface brightness objects have in 
its triggering mechanism. In Paper I (|L6pez-Sanchez fc Estebanll2008l ) we described the motivation of this work, com- 
piled the list of the 20 analysed WR galaxies (Table 1 of Paper I), the majority of them showing several sub-regions or 
objects within or surrounding them ; and presented the results of the optical/ NIR broad-band and Ha photometry In 
Paper II (Lopez-Sanchez & Esteban 2009) we presented the results of the analysis of the intermediate resolution long-slit 
spectroscopy of 16 WR galaxies of ou r samp le - the results for the other four galaxies were pu b lished separately, see 
lLopez-Sanchez. Esteban fc Rodriguez] ([2004al lbT) ; lLopez-Sanchez. Esteban fc Garcfa-Roiasl (|2006l ); lLopez-Sanchez et al.l 
(|2007f) . In many cases, two or more slit positions were used to analyse the m ost interesting zones, knots or mo rphological 
structures belonging to each galaxy or even surrounding objects. Paper III (jLopez-Sanchez fc Estebanl^OlOal) presented 
the analysis of the O and WR stellar populations within these galaxies. Paper IV (jLopez-Sanchez fc Esteban! l2010bT ) 
globally compile and analyse the optical/TVZR photometric data and study the physical and chemical properties of the 
ionized gas within our galaxy sample. The results shown in this paper haven been already published in Paper IV. The final 
paper of the series (Paper V; Lopez-Sanchez 2010) compiles the properties derived with data from other wavelengths 
(UV, FIR, radio, and X-ray) and complete a global analysis of all available multiwavelength data of our WR galaxy 
sample. We have produced the most comprehensive data set of these galaxies so far, involving multiwavelength results 
and analysed according to the same procedures. 



2. Empirical calibrations of the oxygen abundance 

When the spectrum of an extragalactic H II region does not show the [O ill] A4363 emission line or other auroral lines 
that can be used to derive T e , the so-called empirical calibrations are applied to get a rough estimation of its metallicity. 
Empirical calibrations are inspired partly by photo-ionization models and partly by observational trends of line strengths 
with galactocentric distance in gas-rich spirals, which are believed to be due to a radial abundance gradient with abun- 
dances decreasing outwards. In e xtragalactic objects, the usefulness of the empirical methods goes beyond the derivation 
of abundance gradients in spirals ()Pilvugin. Vflchez fc Con tini 2004f) , as these methods find application in chemical a bun- 
dance studies of a variety of objects, including low-surface brightness galaxies (|de Narav. McGaugh fc de Blokl l2004') and 
star-forming galaxies at intermediate and high redshift, where the advent of 8-10 m class telescopes has made it possible 
to extend observations (e.g., Teplitz et al. 2000, Pettini et al. 2001; Kobulnicky et al 2003; Lilly, Carollo & Stockton 
2003; Steidel et al. 2004; Kobulnicky & Kewley 2004; Erb et al. 2006). 

As the brightest metallic lines observed in spectra of H n regions are those involving oxygen, this element has been 
extensively used to get a suitable empirical calibration. Oxygen abundance is important as one of the fundamental 
characteristics of a galaxy: its radial distribution is combined with radial distributions of gas and star surface mass 
densities to constrain models of chemical evolution. Parameters defined in empirical calibrations evolving bright oxygen 
lines are 

J([OIII])A4959 + 7([OIII])A5007 , , 

R3 = H/3 ' (1) 

J([OII])A3727 

R2 = — w — ' {) 

R23 = R3 + R2, (3) 
P = | L , (4) 

y = lo g f 2 =lo gj ^— i . (5) 

iJensen. Strom fc Strom! |l976!) presented the first explorati on in this method considering the R3 index, which con- 
siders the [O III] AA4959,5007 emission lines. However, were iPagel et "ail (|1979[ ) who introduced the most widely used 
abundance indicator, the R23 index, which also included th e bright [Oil] A3727 emi s sion line. Since then, many stud- 
ies have been performed to refine the calibration of R^ (lEdmunds fc Pagell [1985 McCall. Rvbski fc Shields] [l985: 



Dopita fc Evanslll986t iTorres-Peimbert et al.lll989l: lMcGaughlll99lUZaritskv. Kennicutt fc Huchralll994l: iPilvuginlboO^ 



2001allbl : lKewlev fc Dopitall2002l:lKobulnickv fc Kewlevll2004lPilvugin fc Thuanll2005l:lNagao. Maiolino fc Marconill200' 
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Fig. 1. Empirical calibrations of oxygen abundance using the R23 parameter. Note that they are bi-valuated. The dashed zone 
indicates the region with higher uncertainties in O/H. The empirical calibrations plotted in the figure are: EP94: Edmund & Pagel 
(1984); M91: McGaugh (1991) using y=0 (R 2 = R 3 ); P01: Pilyugin (2001) using P = 0.5 (R 2 = R3); (KD02+KK04): Kewley & 
Dopita (2002) using the formulation of Kobulnicky & Kewley (2004) assuming q = 7.5 x 10 7 cm s _1 ; N06: Nagao et al. (2006) 
using their cubic fit to R23 . 



The most successful are the calibrations of lMcGaughl (ll99l|) andlKewlev fc Dopital (l2002|). which are based on photoion- 
ization models, and the empirical relations provided by iPilvuginl (|2001allbf ) and iPilvugin fc Thuanl (pOOl . Both kinds 
of calibrations improve the accuracy by making use of the [O ill] /[On] ratio as ionization parameter, which accounts 
for the large scatter found in the R23 v ersus oxygen abundance calibration, which is larger than observational errors 
(jKobulnickv. Kennicutt fc Pizagno| [l999) . Figure [T] shows the main empirical calibrations that use the -R23 parameter. 

The main problem associated with the use of R23 parameter is that it is bivaluated, i.e., a single value of R23 can be 
caused by two very different oxygen abundances. The reason of this behaviour is that the intensity of oxygen lines does 
not indefinitely increase with metallicity. Thus, there are two branches for each empirical calibration (see Fig. [T]): the 
low-metallicity regime, with 12+log(0/H)<8.1, and the high-metallicity regime, with 12+log(0/H)>8.4. That means 
that a very large fraction of the star-forming regions lie in the ill-defined turning zone around 12+log(O/H)~8.20, where 
regions with the same R23 value have oxygen abundances that differ by almost an order of magnitude. Hence, additional 
information, such as the [Nil] /Ha or the [Oii]/[Oiii] ratios, is needed to break the degeneracy between the high and low 
branches (i.e., Kewley & Dopita, 2002). Besides, the R23 method requires that spectrophotometric data are corrected by 
reddening, which effect is crucial because [0 11] and [O ill] lines have a considerably separation in wavelength. 

Here we list all empirical calibrations that were considered in this work, compiling the equations needed to derive the 
oxygen abundance from bright emission line ratios following every method. 

Edmund & Pagel (1984) : Although th e R23 p arameter was firstly proposed bv lPagel et al.l(|l979h . the first empirical 
calibration was given bv lEdmunds fc Pagell (|1984h . 

12 + log(0/H) up = 8.76 - 0.69 log R 23 , 

12 + \og(0/R) low = 6.43 + 1.67 log R 23 , (6) 

with the limit between the lower and the upper branches at 12+log(O/H)~8.0. 

McCall, Rybski &: Shields (1985) presented an empirical calibration for oxygen abundance using the R23 param- 
eter, only valid for 12+log(0/H)>8.15. However, they did not give an analytic formulae but only listed it numerically 
(see their Table 15). The four-order polynomical fit for their values gives the following relation: 

12 + log(0/H) up = 9.32546 - 0.360465a; + 0.203494a; 2 +0.278702a; 3 - 1.36351x 4 , (7) 
with x = log R23 ■ 

Zaritzky, Kennicutt &: Huchra (1994) provided a simple analytic relation between oxygen abundance and R23: 
12 + log(0/H) up = 9.265 - 0.33a; - 0.202x 2 - 0.207a 3 - 0.333a 4 . (8) 
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Their formula is an average of three previous calibrations: lEdmunds fc Pagell (|l984j) . McCall et al. (1985) and Dopita & 
Evans (1986). Following the authors, this calibration is only suitable for 12+log(O/H)>8.20, but perhaps a more realistic 
lower limit is 8.35. 

McGaugh (1991) calibrated the relationship between t he P 2 3 ratio and g as-phase oxygen abundance using Hn 
region models derived from the photoionization co de Cloudy ((Ferland et al.lll998 ). McGaugh's models include the effects 



of dus t and variations in ionization parameter, v. iKobulnickv et al.l ( 1999(1 1 give analytical expressions for the iMcGaughl 



(1991) calibration based on fits to photoionization models; the middle point between both branches is 12+log(0/H)~8.4 
12 + log(0/H)„ p = 7.056 + 0.767s + 0.602a; 2 - 2/(0.29 + 0.332a; - 0.331a; 2 ), (9) 
12 + log(0/H) ;ou , = 9.061 - 0.2a; - 0.237a; 2 - 0.305a; 3 - 0.0283a; 4 

-?/(0.0047 - 0.0221a; - 0.102a; 2 - 0.0817a; 3 - 0.00717a; 4 ). (10) 

Pilyugin (2000) found that the previous calibrations using the R23 parameter had a systematic error depending on 
the hardness of the ionizing radiation, suggesting that the excitation parameter, P, is a good indicator of it. In several 
papers, Pilyugin performed a detailed analysis of the observation al data combined with photoionization m odels to obtain 
empirical calibrations for the oxygen abundance. iPilvuginl (|2000h confirmed the idea of IMcGaughl (|l99ll ) that the strong 
lines of [0 11] and [O ill] contain the necessary information for the determination of accurate abundances in low-metallicity 
(and may be also in high-metallicity) H 11 regions. He used new observational data to propose a linear fit involving only 
the R23 parameter, 

12 + log(0/H)„ p = 9.50 -1.40 log R 23 , (11) 
12 + log(0/H) io „ = 6.53 + 1.40 log R 23 , (12) 

assuming a limit of 12 +log(O/H)^8.0 between th e two b r anche s. This calibration is close to that given by 
lEdmunds k, Pagel (|1984D ; it has the same slope, but IPilvuginl ([2000) is shifted towards lower abundances by around 
0.07 dex. However, this new relation is not sufficient to explain the wide spread of observational data. Thus, Pilyugin 
(2001a) give the following, more real and complex, calibration involving also the excitation parameter P: 

10 ,, (nm , P 23 + 54.2 + 59.45P + 7.31P 2 

12 + log(0/H) up - 6 01 + 6 71 p + 371 p 2 + 243il . 23 - (13) 

This is the so-called P-method, which can be used in moderately high-metallicity H 11 regions (12+log(Q/H)^ 8.3) . Pilyugin 
used two-zone models of H 11 regions and assumed the T e (On) - T e (Oin) relation from iGarnettl (|l992f ). For the low 
metallicity branch, Pilyugin (2001b) found that 

12 + log(0/H) iou) = 6.35 + 1.45 lo gj R 23 - 1.741ogP. (14) 

Pilyugin estimates that the precision of oxygen abundance determination with this method is around 0.1 dex. 

Pilyugin & Thuan (2005) revisited these calibrations including more spectroscopic measurements of Hn regions 
in spiral and irregular galaxies with a measured intensity of the [O in] A4363 line and recalibrate the relation between 
the oxygen abundance and the P23 and P parameters, yielding to: 

R 23 + 106.4 + 106.8P-3.40P 2 
12 + log(0/HW, = 17-72 + 6 . 60 p + 6 . 95JM _ . 302B23 . (15) 

, n/m P23 + 726.1 + 842.2P + 337.5P 2 

12 + log(0/H)„ p = 85 . 96 + 8276P + 43 . 98p2 + L793 p 23 - (16) 

Kewley & Dopita (2002) used a combination of stellar population synthesis and photoionization models to develop 
a set of ionization parameters and abundance diagnostic based only on the strong optical emission lines. Their optimal 
method uses ratios of [Nil], [On], [Om], [Sn], [Sin] and Balmer lines, which is the full complement of strong nebular 
lines accessible from the ground. T hey also recommend proc edures for the derivation of abundances in cases where only 
a subset of these lines is available. Ke wlev &; Dopita] (|2002t) models start with the assumption that P 2 3, and many of 
the other emission-line abundance diagnostics, also depends on the ionizat ion parameter q = c x U, that has units 
of cm s _1 . They used the stellar population synthesis c odes Starburst 99 (jLeitherer et al.l[l999t IVazauez &; Leithererl 
2005) and Pegase.2 (jFioc fc Rocca-Volmerangelll997tl to generate the ionizing radiation field, assuming burst mod- 
els at zero age with a Salpeter IMF and lower and upper mass limits of 0.1 and 120 M , respectively, with metal- 
licities between 0. 05 and 3 times solar. The ionizing radiation fields were input into the photoionization and shock 
code, Mappings ( Sutherland fc Dopital Il993l ). which includes self-consistent treatment of nebular and dust physics. 



iKewlev &: Dopital (120021) pre viously used these models to simulate the emission-line spectra of H II regions and starburst 

galaxies (jDopita et al.ll2000l) . and are completely described in their s tudy. 

Kobulnicky &: Kewley (2004) gave a parameter ization of the Kewlev fc Dopital (|2002f) P 2 3 method with a form 



similar to that given by IMcGaughl (|1991l ) calibration. IKobulnickv fc Kewlev! ( 20041 ) presented an iterative scheme to 



resolve for both the ionization parameter q and the oxygen abundance using only [Oin], [On] and H/3 lines. The 
parameterization they give for q is 

, /N 32.81 - 1.1532/ 2 + [12 + log(0/H)][- 3.396 - 0.0252/ + 0.1444?/ 2 ] , . 

log: ({7) = — (17) 

6W 4.603 - 0.31192/ - 0.163?/ 2 + [12 + log(0/H)][- 0.48 + 0.0271?/ + 0.02 3 72/ 2 ]' 
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where y = log([0 m]/[0 n]). This equation is only valid for ionization parameters between 5xl0 6 and 1.5 xlO 8 cm s 1 . 
The oxygen abundance is parameterized by 

12 + log(0/H)„ p = 9.72 - 0.777a; - 0.951a: 2 - 0.072a; 3 - 0.811a; 4 - log(g) 

x (0.0737 - 0.0713a; - 0.141a; 2 + 0.0373a; 3 - 0.058a; 4 ), (18) 
12 + log(0/H) /ou , = 9.40 + 4.65a; - 3.17a; 2 - log(g)(0.272 + 0.547a; - 0.513a; 2 ), (19) 

being x = logi?23- The first equation is valid for 12+log(0/H)>8.4, while the second for 12+log(0/H)<8.4. Typically, 
between two and three iterations are required to reach convergence. Following the authors, this parameterization should 
be regarded as an improved , implementation-friendly approach to be preferred over the tabulated R23 coefficients given 
bv lKewlev k Dopital (|2002t) . 

Nagao, Maiolino k Marcani (2006) did not consider any ionization parameter. They merely used data of a large 
sample of galaxies from the Sloan Digital Sky Survey (SDSS; York et al. 2000) to derive a cubic fit to the relation between 
i?23 and the oxygen abundance, 

logi? 23 = 1-2299 - 4.1926y + 1.0246?/ 2 - 0.063169y 3 , (20) 
with y=12+log(0/H). 

Besides i?23, additional parameters have been used to derive metallicities in star-forming galaxies. Without other 
emission lines, the N2 parameter, which is defined by 

(21) 

Ha 

can be used as a crude estimator of metallicity. However, we note that the [Nil] /Ha ratio is particularly sen- 
sitive to shock excitation or a hard radia tion field from an AGN. The A2 parameter was firstly suggested by 
IStorchi-Bergmann. Calzetti k Kinnevl (|l994t ). who gave a tentative calibration of the oxygen abundance using this pa - 
rameter. This ca l ibrati on h as been revisited by Ivan Zee. Salzer k Havnesl (|1998fklDenicol6. Terlevich k Terlevichl (|2002l) ; 
iPettini k Pagell (|2004l ) and lNagao et al.l (|2006h . The Denicolo et al. (2002) calibration is 

12 + log(0/H) = 9.12 + 0.73^2, (22) 

which considerably improves the previous relations because of the inclusion of an extensible sample of nearby extragalactic 
H II regions. The uncertainty of this method is ~0.2 dex because N2 is sensitive to ioniz ation and O/N varia tions, 
so strictly speaking it should be used mainly as an indicator of galaxy-wide abundances. iDenicolo et al.l (|2002|) also 
compared their method with photoionization models, concluding that the observed A^ is consistent with nitrogen being 
a c ombination of both pri mary and secondary origin. 

IPettini k Pagell (|2004f) revisited the relation between the N2 parameter and the oxygen abundance including new 
data for the high- and low-metallicity regimen. They only considered those extragalactic H 11 regions where the oxygen 
values are determined either via the T e method or with detailed photoionization modelling. The linear fit to their data is 

12 + log(0/H) = 8.90 + 0. 57N 2 , (23) 

which has both a lower slope and zero-point that the fit given bv lDenicolo et ahl ((2002). A somewhat better relation is 
provided by a third-order polynomical fit of the form 

12 + log(0/H) = 9.37 + 2.032A 2 + 1.26(Af 2 ) 2 + 0.32(A 2 ) 3 , (24) 

valid in the range —2.5 < N2 < — 0.3. lNagao et ahl (|2006F) also provided a relation between N2 and the oxygen abundance, 
their cubic fit to their SDSS data yields 

log N 2 = 96.641 - 39.941?/ + 5.2227?/ 2 - 0.22040?/ 3 , (25) 
with ?/=12+log(0/H). 

Pettini k Pagel (2004) revived the O3N2 parameter, previously introduced bv lAlloin et al.l (|1979|) and defined by 
^ Ar , [O III] A5007/H£ 

Q37V2 " l0S [NII]A6583/Ha - (26) 
IPettini k Pagell (|2004l ) derived the following least-square linear fit to their data: 

12 + log(0/H) = 8.73 - O.32O3A2. (27) 

INagao et al.l (|2006l ) also revisited this calibration and derived a cubic fit between the O3N2 parameter and the oxygen 
abundance, 

log 3 A 2 = -232.18 + 84.423y - 9.9330?/ 2 + 0.37941?/ 3 , (28) 
with y=12+log(0/H). 
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Table 1. List of the parameters used to compute the oxygen abundance in all regions with a direct estimation of T e using empirical 
calibrations. 



Region 


1 1 

^23 


f = K3/K23 


y =iog(«3/K2) 


iV2 




QKD02o a 


HCG 31 AC 


5.42 


0.571 


0.125 


0.104 


1.349 


3.76E+07 


HCG 31 B 


7.93 


0.408 


-0.162 


0.101 


1.381 


4.91E+07 


HCG 31 E 


7.12 


0.511 


0.020 


0.090 


1.486 


7.40E+07 


HCG 31 Fl 


8.91 


0.819 


0.656 


0.034 


2.201 


5.78E+07 


HCG 31 F2 


7.60 


0.724 


0.418 


0.036 


2.064 


6.28E+07 


HCG 31 G 


8.20 


0.499 


-0.002 


0.106 


1.462 


6.96E+07 


Mkn 1199 C 


7.69 


0.809 


0.627 


0.131 


1.555 


1.55E+08 


Haro 15 A 


9.73 


0.884 


0.881 


0.027 


2.378 


8.55E+07 


Mkn 5 Al 


7.58 


0.748 


0.473 


0.051 


1.915 


6.96E+07 


Mkn 5 A2 


8.19 


0.702 


0.372 


0.049 


1.944 


1.72E+08 


IRAS 08208+2816 C 


7.77 


0.793 


0.583 


0.129 


1.558 


8.55E+07 


POX 4 


10.68 


0.906 


0.986 


0.015 


2.697 


1.05E+08 


UM 420 


6.45 


0.649 


0.268 


0.099 


1.497 


4.81E+07 


SBS 0926+606A 


7.40 


0.811 


0.632 


0.026 


2.227 


6.68E+07 


SBS 0948+532 


8.85 


0.874 


0.843 


0.022 


2.430 


2.54E+08 


SBS 1054+365 


9.33 


0.893 


0.920 


0.020 


2.503 


9.10E+07 


SBS 1211+540 


7.22 


0.892 


0.918 


0.008 


2.788 


1.16E+08 


SBS 1319+579A 


9.92 


0.908 


0.996 


0.014 


2.671 


1.05E+08 


SBS 1319+579B 


7.13 


0.722 


0.415 


0.046 


1.922 


6.15E+07 


SBS 1319+579C 


7.11 


0.710 


0.389 


0.052 


1.860 


5.91E+07 


SBS 1415+437C 


5.22 


0.783 


0.558 


0.015 


2.301 


5.91E+07 


SBS 1415+437A 


4.86 


0.810 


0.629 


0.012 


2.370 


5.44E+07 


III Zw 107 A 


7.13 


0.701 


0.369 


0.100 


1.573 


5.78E+07 


Tol 9 INT 


4.58 


0.689 


0.345 


0.252 


0.973 


4.16E+07 


Tol 9 NOT 


4.78 


0.629 


0.230 


0.287 


0.894 


3.39E+07 


Tol 1457-262A 


9.89 


0.773 


0.532 


0.033 


2.236 


9.91E+07 


Tol 1457-262B 


9.00 


0.792 


0.582 


0.020 


2.417 


1.41E+08 


Tol 1457-262C 


8.88 


0.669 


0.359 


0.036 


2.099 


7.16E+07 


ESO 566-8 


5.17 


0.505 


0.008 


0.414 


0.693 


3.19E+07 


NGC 5253 A 


9.20 


0.851 


0.756 


0.102 


1.754 


6.82E+07 


NGC 5253 B 


9.38 


0.856 


0.775 


0.086 


1.841 


7.11E+07 


NGC 5253 C 


8.03 


0.773 


0.532 


0.041 


2.046 


2.60E+08 


NGC 5253 D 


7.67 


0.527 


0.048 


0.079 


1.582 


7.72E+07 



a Value derived for the q parameter (in units of cm s 1 ) obtained using the optimal calibration given bv lKewlev fc Dopital [12002, 1 . 



Other important empirica l calibrations that were n ot used in this study involve the 6*23 parameter, intro- 
duced bv IVilchez fc Estebanl (|1996t ) and revisited by iDiaz fe Perez-Monterol (|2000T) : lOev fe Shields! (j2000l ) and 
iPerez-Montero fe Diazl ([20051 ). In the last years, b right emission line ratios such as [Ar iil]/[0 ill] and [S III]/ [O III] 
(jStasihskal I2006T ) or~[Ne Hi] / [O ill] and [O III] /[Oil] (iNaeao et al.l I2006T) ha ve been explored as indicators of the oxy- 
gen abundance in Hli regions and starburst galaxies. iPeimbert et al. ( 20071) suggested to use the oxygen recombination 
lines to get a more precise estimation of the oxygen abundance. Nowadays, there is still a lot of observational and the- 
oretical work to do involving empirical calibrations (see recent review by Kewley & Ellison 2008), but these methods 
should be used o nly for objects whose H II regions have the same structural properties as those of the calibrating samples 
(|Stasihsldl2009h . 



3. Comparison with empirical calibrations 

We used the dat a of the 31 regions for whic h we ha ve a direct estimate of T e and, hence, a direct estimate of the oxygen 
abundance -see lLopez-Sanchez fe Estebanl (|2010bf) and their Table 3- , to check the reliability of several empirical 
ca librations. A recent revie w of 10 metallicity calibrations, including theoretical and empirical methods, was presented 
bv lKewlev fc Ellisonl ((2008), but previous section gives an overview of the most common empirical calibrations and defines 
the typical parameters that are used to estimate the oxygen abundance following these relations. These parameters are 
ratios between bright emission lines, the most commonly used are R23, P, y, N2, and O3-/V2 T able [J lists the values of all 
these parameters derived for each region with a direct estimate of the oxygen abundance see Lopez-Sanchez & Esteban 
(2009) for details-. Table [T] also in cludes the value d e rived for the q parameter (in units of cm s _1 ) obtained from 
the optimal calibration provided bv iKewlev fc D opita (2002). The results for the oxygen abundances derived for each 
object and empirical calibration are listed in Table [2] This table also indicates the branch (high or low metallicity) 
considered in each region when using the R23 parameter although, as is clearly specified in the table, for some objects 
with 8.00<12+log(O/H)<8.3 we assumed the average value found for the lower and upper branches. 

Looking at the data compiled in Table [21 the huge ra nge of oxygen abundance found for the same object using 
different calibrations is evident. As lKewlev fc Ellisonl (|2008f) concluded, it is critical to use the same metallicity calibration 
when comparing properties from different data sets or investigate luminosity-metallicity or mass-metallicity relations. 
Furthermore, abundances derived with such strong-line methods may be significantly biased if the objects under study 
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Table 2. Results of the oxygen abundance, in the form 12+log(0/H), for objects with a direct estimation of the metallicity, 
considering several empirical calibrations. 



JrtGglOU 


ISranch 


T 




EP84 MRS85 


1\/TQ1 


ZKH94 POO 


P01 a 


PT05 6 


KD02 


KK04 D02 PP04a PP04b PP04c N06a c 


N06b N06c 










R23 


R23 


R23, y 


R23 


R23 R23, PR23, PR23, ?R23, q N 2 


N 2 


N 2 


O3N2 


R23 


N 2 


O3N2 


HCG 31 AC 


H 


8.22+0 


05 


8.25 


8.89 


8.67 


8.74 


8.47 


8.15 


8.09 


7.99 


8.12 


8.40 


8.34 


8.29 


8.30 


8.05 


8.16 


8.22 


HCG 31 B 


H 


8.14+0 


08 


8.14 


8.48 


8.29 


8.44 


8.24 


8.22 


8.12 


8.41 


8.44 


8.39 


8.33 


8.28 


8.29 


8.07 


8.16 


8.20 


HCG 31 E 


H 


8.13+0 


09 


8.17 


8.62 


8.14 


8.53 


8.31 


8.18 


8.13 


8.19 


8.32 


8.35 


8.30 


8.26 


8.25 


8.07 


8.11 


8.15 


HCG 31 Fl 


L 


8.07+0 


06 


8.02 


8.30 


8.13 




7.86 


8.12 


7.99 


8.46 


8.33 


8.05 


8.07 


8.09 


8.03 




7.81 


7.67 


HCG 31 F2 


L 


8.03+0 


10 


7.90 


8.54 


8.06 


8.48 


7.76 


8.13 


7.95 


8.19 


8.27 


8.06 


8.07 


8.10 


8.07 


8.07 


7.83 


7.79 


HCG 31 G 


H 


8.15+0 


08 


8.13 


8.43 


8.26 


8.40 


8.22 


8.11 


8.17 


8.31 


8.42 


8.41 


8.34 


8.29 


8.26 


8.07 


8.17 


8.16 


Mkn 1199 C 


H 


8.75+0 


12 


9.37 


9.26 


9.00 


9.18 


9.19 


8.71 


8.54 


9.14 


9.14 


8.92 


8.74 


8.90 


8.81 


9.18 


8.78 


8.94 


Haro 15 A 


H 


8.10+0 


06 


8.08 




8.14 




7.91 


8.12 


8.12 


8.48 


8.34 


7.98 


8.01 


8.05 


7.97 




7.74 


7.38 


Mkn 5 Al 


L 


8.07+0 


07 


7.90 


8.54 


8.04 


8.48 


7.76 


8.13 


8.13 


8.19 


8.26 


8.18 


8.17 


8.16 


8.12 


8.07 


7.94 


7.89 


Mkn 5 A2 


L 


8.08+0 


07 


7.95 


8.43 


8.14 


8.41 


7.81 


7.92 


8.17 


8.18 


8.33 


8.16 


8.15 


8.15 


8.11 


8.07 


7.92 


7.87 


IRAS 08208+2816 


H 


8.33+0 


08 


8.15 


8.50 


8.55 


8.46 


8.25 


8.42 


8.35 


8.35 


8.25 


8.47 


8.39 


8.34 


8.23 


8.35 


8.23 


8.11 


POX 4 


L 


8.03+0 


04 


8.15 




8.20 




7.97 


7.92 


8.06 


8.48 


8.40 


7.78 


7.86 


7.91 


7.87 




7.53 




UM 420 


L 


7.95+0 


05 


7.78 


8.73 


7.98 


8.61 


7.66 


7.85 


7.86 


8.02 


8.16 


8.20 


8.39 


8.33 


8.28 


7.57 


8.15 


8.14 


SBS 0926+606A 


L 


7.94+0 


08 


7.88 


8.57 


7.97 


8.50 


7.75 


7.77 


7.80 


8.17 


8.20 


7.97 


8.00 


8.05 


8.02 


7.71 


7.73 


7.64 


SBS 0948+532 


L 


8.03+0 


05 


8.01 


8.31 


8.06 




7.86 


7.82 


8.10 


8.34 


8.28 


7.91 


7.95 


8.01 


7.95 




8.01 




SBS 1054+365 


L 


8.00+0 


07 


8.05 


8.21 


8.09 




7.89 


7.84 


7.91 


8.48 


8.30 


7.87 


7.93 


7.98 


7.93 




7.63 




SBS 1211+540 


L 


7.65+0 


04 


7.86 


8.60 


7.85 


8.52 


7.73 


7.68 


7.65 


8.02 


8.10 


7.58 


7.70 


7.69 


7.84 


7.68 


7.31 




SBS 1319+579A 


L 


8.05+0 


06 


8.09 




8.13 




7.93 


8.11 


8.11 


8.48 


8.33 


7.77 


7.85 


7.90 


7.88 




7.52 




SBS 1319+579B 


L 


8.12+0 


10 


7.85 


8.62 


8.01 


8.53 


7.72 


8.13 


8.12 


8.13 


8.23 


8.14 


8.14 


8.14 


8.11 


8.07 


7.90 


7.89 


SBS 1319+579 C 


L 


8.15+0 


07 


7.85 


8.62 


8.02 


8.53 


7.72 


8.13 


8.13 


8.12 


8.23 


8.18 


8.17 


8.16 


8.13 


8.06 


7.94 


7.93 


SBS 1415+437 C 


L 


7.58+0 


05 


7.63 


8.91 


7.72 


8.76 


7.53 


7.57 


7.55 


7.86 


7.99 


7.79 


7.86 


7.92 


7.99 


7.39 


7.55 


7.55 


SBS 1415+437 A 


L 


7.61+0 


06 


7.58 


8.96 


7.64 


8.80 


7.49 


7.50 


7.48 


7.82 


7.92 


7.72 


7.81 


7.86 


7.97 


7.34 


7.48 


7.41 


III Zw 107 


H 


8.23+0 


09 


8.17 


8.62 


8.57 


8.53 


8.31 


8.40 


8.35 


8.13 


8.24 


8.39 


8.33 


8.28 


8.23 


8.46 


8.15 


8.10 


Tol 9 INT 


H 


8.58+0 


15 


8.30 


9.00 


8.76 


8.84 


8.58 


8.61 


8.55 


8.95 


8.90 


8.68 


8.56 


8.54 


8.42 


8.77 


8.46 


8.40 


Tol 9 NOT 


H 


8.55+0 


16 


8.29 


8.97 


8.73 


8.81 


8.55 


8.56 


8.50 


8.94 


8.88 


8.72 


8.59 


8.59 


8.44 


8.75 


8.51 


8.44 


Tol 1457-262A 


L 


8.05+0 


07 


8.09 




8.26 




7.92 


8.11 


8.20 


8.58 


8.42 


8.04 


8.06 


8.09 


8.02 




7.80 


7.63 


Tol 1457-262B 


L 


7.88+0 


07 


8.55 




8.21 




7.87 


7.91 


8.21 


8.38 


8.29 


7.88 


7.93 


7.99 


7.96 




7.60 




Tol 1457-262C 


L 


8.06+0 


11 


8.48 




8.21 




7.88 


8.00 


8.24 


8.48 


8.37 


8.07 


8.08 


8.10 


8.06 




7.83 


7.77 


ESO 566-8 


H 


8.46+0 


11 


8.27 


8.92 


8.68 


8.77 


8.50 


8.44 


8.38 


8.92 


8.84 


8.84 


8.68 


8.76 


8.51 


8.70 


8.66 


8.53 


NGC 5253 A 


H 


8.18+0 


04 


8.09 


8.24 


8.13 




8.15 


8.11 


8.13 


8.53 


8.33 


8.40 


8.34 


8.28 


8.17 




8.16 


8.00 


NGC 5253 B 


H 


8.19+0 


04 


8.09 


8.21 


8.14 




8.14 


8.11 


8.13 


8.48 


8.34 


8.34 


8.29 


8.25 


8.14 




8.10 


7.94 


NGC 5253 C 


L 


8.28+0 


04 


8.14 


8.46 


8.53 


8.42 


8.23 


8.38 


8.32 


8.67 


8.63 


8.11 


8.11 


8.13 


8.08 


8.30 


7.87 


7.80 


NGC 5253 D 


L 


8.31+0 


07 


8.15 


8.52 


8.19 


8.47 


8.26 


8.23 


8.17 


8.32 


8.37 


8.31 


8.27 


8.23 


8.22 


8.37 


8.07 


8.10 



NOTE: The empirical calibrations and the parameters used for each of them are: EP84: Edmunds & Pagel (1984) that involves the R23 
parameter; MRS85: McCall, Rybski & Shields (1985) using R 23 ; M91: McGaugh (1991) using R 23 and y; ZKH94: Zaritzky, Kennicutt & 
Huchra (1994) using i? 2 3i P00: Pilyugin (2000) using R23; P01: Pilyugin (2001a,b) using R23 and P; KD02: Kewley & Dopita (2002) using 
R23 and q; KK04: Kobulnicky & Kewley (2004) using R 23 & <?; D02: Denicolo, Terlevich & Terlevich (2002) using the N 2 parameter; PP04: 
Pettini & Pagel (2004), using (a) 7V2 with a linear fit, (b) N 2 with a cubic fit, (c) the O3N2 parameter; N06: Nagao et al. (2006) using the 
cubic relations involving the R23 (a), N 2 (b) and O3./V2 (c) parameters. The value compiled in the column labeled T e is the oxygen 
abundance derived by the direct method. 

a The value listed for P01 is the average value between the high- and the low-metallicity branches for objects with 7.90<12+log(O/H)<8.20. 
b The value listed for PT05 is the average value between the high- and the low-metallicity branches for objects with 
8.05<12+log(O/H)<8.20. 

c The value listed for N06 is the average value between the high- and the low-metallicity branches for objects with 8.00<12+log(O/H)<8.15. 



have different structural p roperties (hardn ess of the ionizing radiation field, morphology of the nebulae) than those used 
to calibrate the methods (Stasihska 200jl). 

Figures [2] and [3] plots the ten most common calibrations and their comparison with the oxygen abundance obtained 
using the direct method. We performed a simple statistic analysis of the results to quantify the goodness of these empirical 
calibrations. Table [3] compiles the average value and the dispersion (in absolute values) of the difference between the 
abundance given by empirical calibration and that obtained using the direc t method. We check that the empirical 
calibration that provides the best results is that proposed by iPilvuginl (|2001al fbT) , which gives oxygen abundances very 
close to the direct values (the differences are lower than 0.1 dex in the majority of the objects), and furthermore it 
possesses a low dispersion. We note however that the largest d ivergences found u s ing th is calibration are in the low- 
metallicity regime. The update of this calibration presented by iPilvuein fe Thuanl (|2005t) seems to partially solve this 
problem, the abundances provi ded by this calibration also agree very well with those derived following the direct method. 
We therefore conclude that the IPilvuein fe Thuanl (|2005f) calibration is nowadays the best suitable method to derive the 
oxygen abundance of star-forming galaxies when auroral lines are not observed. 

On the other hand, the results given by the empirical calibrations provided by iMcGaughl ()199lD . iKewlev fc Dopital 
(|2002t) and lKobulnickv fc Kewlevl (|2004l) . that are based on photoionization models, are systematically higher than the 
values derived from the direct method. This effect is even more marked in the last two calibrations, which usually are 
between 0.2 and 0.3 de x higher than th e exp ected values. These empir ical calibrati ons also h a ve a higher dispersion 
than that estimated for IPilvuginl (|2001allbT ) or iPilvugin fe Thuanl (|2005l) calibrations. lYin et all (|2007l) also found high 
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12+log (0/H)-Te 



7.8 8.0 8.2 8.4 

12+log (0/H)-Te 



Fig. 2. Comparison between the oxygen abundances derived using the direct method (T e , always plotted in the x — axis) with th ose 
estima ted using six different emp irical ca librations that consider the R 23 param eter: M91:lMcGaughl (|1991T> : P01 :lPilvuginl (l200f dlbf) ; 
KD02: iKewlev fc Dopital (|2002l ): KK04: iKobulnickv fc Kewlevl ^OO^ : PT05: iPilvuein fc Thuanl l|2005l ): N06: iNaeao et al.l l|2006l i. 
The bottom panel of each diagram indicates the difference between empirical and direct data. 



discrepancies when comparing the theo retical m e tallicit ies usi ng the theoretical models of iTremonti et af] (|2004f ) with 
the r o -based metallicites obtained from iPilvugml t)2001allbh and lPilvugin fc Thuanl (|2005l) . 

One of the possible explanations for the different metallicities obtained between the direct method and those de- 
rived from the empirical calibrations based on photoionization models are temperature fluctuations in the ionized gas. 
Temperature gradients or fluctuations indeed cause the true metallicities based on the T c -method to be underestimated 
(i.e. Peimbert 1967; Stasinska 2002,2005; Peimbert et al. 2007). Temperature fluctuations can also explain our results 
for NGC 5253 (jLopez-Sanchez et alJ[2007l) : the ionic abundances of ++ /H + and C ++ /H + derived from recombination 
lines are systematically 0.2 - 0.3 dex higher than those determined from the direct method - based on the intensity ra- 
tios o f collisionally excited line s. This abundance d iscrepancy has been also found in Galactic (|Garcfa-Roias fc Estebanl 
2007) and other extragalactic i Esteban et alJ[2009fl H 11 regions and interestingly this discrepancy is in all cases of the 
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'. . *. fbjpj 









12+log (0/H)-Te 



7.8 8.0 8.2 8.4 

12+log (0/H)-Te 




7.8 8.0 8.2 8.4 

12+log (0/H)-Te 



Fig. 3. Comparison between the oxygen abunda nces derived using the direct me t hod ( T e , in the x—axis) with those com puted using 
the calibrations using the N 2 parameter -D02 : iDenicolo. Terlevich fc Terlevichl (|2002r ): PP04b: iPettini fc Pagell (12004T) considerin g 
their cubic fit to N2 ; N06: iNagao et al.l (|2006T )- and the O3N2 parameter following the cubic fit provided bv lNagao 
The bottom panel of each diagram indicates the difference between empirical and direct data. 



same order as the differences between abundances derived from the direct methods and empirical calibrations based on 
photoionization models. 

The conclusion that temperature fluctuations do exist in the ionized gas of starburst galaxies is very important for the 
analysis of the chemical evolution of galaxies and the Universe. Indeed, if that is correct, the majority of the abundance 
determinations in extragalactic objects following the direct method, including those provided in this work, have been 
underestimated by at least 0.2 to 0.3 dex. Deeper observations of a large sample of star- forming galaxies -that allow 
us to detect the faint recombination lines, such as those provided by lEsteban et al.l (|2009f) - and mor e theoretical work 
-inclu ding a better understanding of the photoionization models, such as the analysis provided by iKewlev fe Ellison] 
(2008) are needed to confirm this puzzling result. 

On the other hand, we have checked the validity of the recent relation provided by Nae ao et al.l ((2006) , which merely 
considers a cubic fit between the R23 parameter and the oxygen abundance. This calibration was obtained combining 
data from several large galaxy samples, the m ajority from the SP SS, which includes all kinds of star-forming objects. As 
it is clearly seen in Tableland in Fig. [31 the iNaeao et al.l (|2006f) relation is not suitable to derive a proper estimate of 
the oxygen abundance for the majority of the objects in our galaxy sample. In general, this calibration provides lower 
oxygen abundances in low-metallicity regions and higher oxygen abundances in high-metallicity regions. Objects located 
in the metallicity range 8.00<12-r-log(O/H)<8.15 have systematically 12-r-log(O/H)N06 ^8.07 because we have to use an 
avera ge value betw e en th e low and the high branches. Furthermore, many of the regions do not have a formal solution 
to the Nagao et al.l (|2006l ) equation, such as the maxim um value for i?23 is 8.39 at 12-|-log(O/ H)=8.07. We consider tha t 
the use of an ionization parameter -P as introduced by iPilvueinl (|2001al fbT) or q as followed bv IKewlev fc Dopital ([2002) 
is fundamental to obtain a real estimate of the oxygen abundance in star- forming galaxies, e specially in o b jects showing 
strong starbursts. In the same sense, the direct method and not the formulae provided by llzotov et al.l ((2006| ) (which 
assumes a low-density approximation in order not to have to solve the statistical equilibrium equations of the +2 ion) 
provides a good approximation to the actual oxygen abundance when the auroral line [O III] A4653 is observed. 

Empirical calibrations considering a linear fit to the N2 ratio (jDenicolo. Terlevich fc Ter levich 20021: IPettini fc Pagell 

12001 give results that are systematically ^0.15 dex higher that the oxygen abundances derived from the direct method. 
The difference is higher at higher metallicities. We do not consider that this trend is a consequence of comparing different 
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Table 3. Results of the comparison between the oxygen abundance given by several empirical calibrations and the oxygen 
abundance derived here following the direct (T c ) method. 



Parameter 






R23 








N 2 




3 N 2 


Calibration" 


P01 


PT05 


N06 M91 


KD02 


KK04 


D02 


PP04 


N06 


PP04 N06 


Average 6 
Notes'* 


0.07 
0.05 
B/A (1) 


0.08 
0.07 

B/A 


0.14 0.15 
0.12 0.11 
(2) S.H. 


0.28 
0.18 
S.H. 


0.27 
0.13 
S.H. 


0.14 
0.10 
S.H. (3) 


0.12 
0.10 

B/A (4) 


0.18 
0.14 
S.L.(5) 


0.12 0.21 
0.10 0.16 
B/A (6) S.L. 



a The names of the calibrations are the same as in Table [2] 

h Average value (in absolute values) of the difference between the abundance given by empirical calibrations and that obtained 
using the direct method. The names of the calibrations are the same as in Table [2] 

c Dispersion (in absolute values) of the difference between the abundance given by empirical calibrations and that obtained using 
the direct method. 

d We indicate if the empirical calibration gives results both below and above the direct value (B/A), if they are systematically 
higher than the direct value (S.H.) or if they are systematically lower than the direct value (S.L.). Some additional notes are: 

(1) Higher deviation in the low branch. 

(2) This calibration provides lower oxygen abundances in low-metallicity regions and higher oxygen abundances in 
high-metallicity regions. 

(3) Systematically higher only for 12+log(0/H)>8.2. 

(4) Higher deviation for 12+log(O/H)<8.0. Considering 12+log(O/H)>8.0, we get average=0.08 and <r=0.06. 

(5) Higher deviation at lower oxygen abundances. 

(6) Higher deviation for 12+log(O/H)<8.0. Assuming 12+log(O/H)>8.0, average=0.09 and a=0.06. 



objects: both IDenicolo et al.l (|2002T ) and lPettini fc Pagell (|2004f ) calibrations are ob tained using a sample of star- forming 
galaxies similar to those analysed in this work, many of which are WR galaxies. IDenicolo et al.l (|2002f) compared the 
N2 ratio with the ionization parameter together with the results of photoionization models and concluded that most of 
th e observed trend of N 2 with the oxygen abundance is caused by metallicity changes. The cubic fit to N 2 performed 
by iPettini fc Pagell (|2004h better reproduces the oxygen abundance, especially in the intermediate- and high-metallicity 
regime (12+ l og(O/ H)>8.0), where it has an average error of ^0.08 dex. However, the cubic fit to -/V 2 provided by 
iNagao et al.l (|2006l) gives systematically lower values for the oxygen abundance than those derived using the direct 
method, having an average error of ~0.18 dex. 

Th e empirical calibration between the oxygen abundance and the O3N2 parameter proposed by IPettini fe Pagell 
(|2004[) gives accepta ble results fo r objec ts with 12+log(O/H)>8.0, with the average error ~0.1 dex. However, the new 
relation provided bv lNagao et al.l ((2006) involving the O3 N2 parameter gives systematically lower values for the oxygen 
abundance. As we commented before, we consider that the lNaeao et al.l (|2006f) calibrations are not suitable for studying 
galaxies with strong star-formation bursts. Their procedures must be taken with caution, galaxies with different ionization 
parameters, different chemical evolution histories, and different star formation histories should have different relations 
between the bright emission lines and the oxygen abundance. This issue is even more important when estimating the 
metallicities of intermediate- and high-redshift galaxies, because the majority of their properties are highly unknown. 

4. Conclusions 

We compared the abundances provided by the direct method with those obtaine d using the most common empiric al 
calibrations in our sample of star-forming regions within Wolf-Rayet galaxies -see lLopez-Sanchez fc Estebanl (|2010bf )-. 
The main conclusions are: 

— The Pilyugin-method of iPilvueinl (2001a. b), which considers the R23 and the P parameters and is updated by 
iPilvugin fc Thuanl (|2005j ). is nowadays the best suit able empirica l calib ration to derive the oxygen abundance of star- 
forming galaxies. The cubic fit to R23 provided by INagao et al.l ((2006) is not valid for analysing these star-forming 
galaxies. 

— The relations between the oxygen abundance and the N2 or the O3N2 parameters provided by IPettini fc Pagel (|2004l) 
give acceptable results for objects with 12+log(O/H)>8.0. 

— The r esults provided by empirical calibrations based on photoionization models (|McGaughlll99ll : iKewlev fc Dopital 
2002; iKobulnickv fc Kewlevl I20&41) are systematically 0.2 - 0.3 dex higher than the values derived from the direct 
method. These differences are of the same order as the abundance discrepancy found between abundances determined 
from recombination and collisionally excited lines of heavy-element ions. This may suggest temperature fluctuations 
in the ionized gas, as they exist in Galactic and other extragalactic H 11 regions. 
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